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We present a framework for relativistic quantum metrology that is useful for both Earth-based and 
space-based technologies. Quantum metrology has been so far successfully applied to design precision 
instruments such as clocks and sensors which outperform classical devices by exploiting quantum 
properties. There are advanced plans to implement these and other quantum technologies in space, for 
instance Space- QUEST and Space Optical Clock projects intend to implement quantum communications 
and quantum clocks at regimes where relativity starts to kick in. However, typical setups do not take into 
account the effects of relativity on quantum properties. To include and exploit these effects, we introduce 
techniques for the application of metrology to quantum field theory. Quantum field theory properly 
incorporates quantum theory and relativity, in particular, at regimes where space-based experiments take 
place. This framework allows for high precision estimation of parameters that appear in quantum field 
theory including proper times and accelerations. Indeed, the techniques can be applied to develop a novel 
generation of relativistic quantum technologies for gravimeters, clocks and sensors. As an example, we 
present a high precision device which in principle improves the state-of-the-art in quantum accelerometers 
by exploiting relativistic effects. 

Quantum technologies are widely expected to bring about many key technological advances this century. 
Experiments in quantum communication are rapidly progressing from table-top to space-based setups. 
For instance, in 2012 a teleportation protocol was successfully performed across a distance of 143 km by 
the group led by A. Zeilinger\ Partly motivated by this success, major space agencies, e.g., in Europe and Canada, 
have invested resources for the implementation of space-based quantum technologies^"^. There are advanced 
plans to use satellites to distribute entanglement for quantum cryptography and teleportation (e.g., the Space- 
QUEST project^) and to install quantum clocks in space (e.g., the Space Optical Clock project^). However, at these 
scales relativistic effects become observable. General relativity provides an effective description of the Universe at 
large length scales — observable effects can thus be expected at the regimes where satellites operate. For instance, 
the Global Positioning System (GPS), a system of satellites used for time dissemination and navigation, requires 
relativistic corrections to determine time and positions accurately^. Cutting-edge experiments are reaching 
relativistic regimes, yet the effects of gravity and motion on quantum technologies are largely unknown. 

The inability to unify quantum theory and general relativity remains one of the biggest unsolved problems in 
physics today. Understanding general relativity at small length scales where quantum effects become relevant is a 
highly non-trivial endeavour that suffers from a scarcity of experimental guidance. Therefore, an alternative 
approach is to study quantum effects at large scales where experiments promise to be within reach in the near 
future^'^. However, in order to exploit quantum properties to measure position and time at scales where relativity 
becomes relevant, it is necessary to work within quantum field theory which provides a description of quantum 
fields in curved space-time^°. It is a semiclassical description in the sense that matter and radiation are quantized 
but the spacetime is classical. As a first step in this direction, it was shown that quantum metrology techniques can 
be applied to measure the Unruh effect at accelerations that are within experimental reach^^"^^ and to estimate the 
curvature of space- time It was also shown that entanglement can be used to determine space- time parameters 
such as the expansion rate of the universe^^. An appropriate quantum field theory approach that includes the 
effects of quantum theory and relativity will enable the development of a new generation of quantum technologies 
for space. Indeed, previous work on relativistic quantum information has already addressed fundamental ques- 
tions about entanglement in quantum field theory where results show that gravity, motion and space-time 
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dynamics can create and degrade entanglement^^. Recent work^^ 
shows that this impHes that relativistic motion produces observable 
effects on quantum communication. These preliminary results are of 
great importance for the space-based quantum experiments that will 
test quantum theory at large scales. 

In this paper we develop a new framework for relativistic quantum 
metrology by incorporating and exploiting relativistic effects in 
quantum parameter estimation. The framework provides the neces- 
sary methodology that will enable the design and production of new 
measurement instruments. As an example we present an accel- 
erometer (see Fig. 1) with a precision that, in principle, improves 
the state-of-the-art in cutting-edge experiments to detect gravity 
anomalies in outer space^^'^^. Our device is based on the fact that 
acceleration produces observable relativistic effects on Bose-Einstein 
Condensates (BEC)^°. The motion-induced transformation of the 
state of the relativistic phononic field on the EEC allows an extremely 
accurate estimation of the magnitude of the acceleration. We use 
quantum metrology tools, in particular, quantum Fisher informa- 
tion^^ to compute optimal bounds to the error of the estimation 
procedure. 

Methods 

Quantum field theory and the covariance matrix formalism. We are interested in 
applying metrology techniques to determine with high precision parameters that 
appear in quantum field theory, for instance accelerations, gravitational field 
strengths, and proper times. In order to do so, we begin by reviewing basic concepts 
from quantum field theory and the covariance matrix formalism. We consider a real, 
scalar quantum field that obeys the Klein- Gordon equation in curved spacetime. It is 
convenient to expand the field in a discrete set of modes = 1, 2, 3, . ..Y^ so that 
we can represent the field and its transformations in the covariance matrix formalism 
which is commonly employed in continuous variable quantum information and 
quantum metrology. We refer the readers to^^'^^ for further details. The functions 
are solutions to the field equation and form a complete set of orthonormal modes with 
respect to a chosen inner product^" denoted by (.,.). The creation and annihilation 
operators associated to the field modes satisfy the commutation relations 
[am,cin] = = 0 [(im,(il] = The vacuum state is defined as the state that 

is annihilated by the operators a„ for all n, i.e., a„|0> = 0. A coordinate transformation 
between different observers, for example, between inertial and accelerated 
observers^^, results in a Bogoliubov transformation between and mode solutions in 
the new coordinate system denoted by The operators associated with 0„ are a„. 

The most general linear transformation between the field operators and a„ is 
given by. 



(1) 



(a) 



(b) 




Figure 1 | (a) General cavity framework: the initial state of two modes of a 
quantum field inside a cavity, represented initially by the covariance matrix 
Gkk' undergoes a relativistic transformation which depends on some 
parameter 0. The transformed covariance matrix dkk' depends on the 
parameter 0, which can be estimated using quantum metrology tools, 
(b) Example: measurement of the acceleration in a BEG setup. 



where = j and = — \J)n->^mj Bogoliubov coefficients. The 

transformed vacuum |6) is annihilated by the new annihilation operators a„ for all n. 
Note that |o) is annihilated by the initial field operators a„ only if all coefficients 
are zero. Indeed particle production manifests when 0. This occurs, for 

instance, in the Unruh effect where the inertial vacuum state is seen as a thermal state 
by uniformly accelerated observers^". Another example of interest is that of a cavity in 
non-inertial motion. The vacuum state of an inertial cavity becomes populated by 
particles after the cavity undergoes non-uniformly accelerated motion^*^. 

Let us now consider the covariance matrix formalism, which has been very useful to 
investigate entanglement in quantum field theory^'^'^^'^*'. In this phase space form- 
alism, for Gaussian states of a bosonic field, all the relevant information about the 
state is encoded in the first and second moments of the field. In particular, the second 
moments are described by the covariance matrix ay = {XjXj + XjXi) — 2{Xi) (Xj), 
where (.) denotes the expectation value and the quadrature operators X, are the 
generalized position and momentum operators of the field modes. In this paper we 
follow the conventions used in^^-^"^, i.e., the operators for the n-th mode are given by 
1 



X: 



. The covariance matrix formalism 



V2 ^ ' 
enables elegant and simplified calculations and has been proven useful to define and 
analyze computable measures of bipartite and multipartite entanglement for 
Gaussian states^^'^*^. 

Every unitary transformation in Hilbert space that is generated by a quadratic 
Hamiltonian can be represented as a symplectic matrix S in phase space. These 
transformations form the real symplectic group Sp(2/i,[R), the group of real i^n X In) 
matrices that leave the symplectic form Q invariant, i.e., SQS^ = Q, where 

Q = 0"^ ^Q, and = ^ J ^ . The time evolution of the field, as well as the 

Bogoliubov transformations, can be encoded in this structure. The symplectic matrix 
corresponding to the Bogoliubov transformation in Eq. (1) can be written in terms of 
the Bogoliubov coefficients as 



\ 

where the hAmn are the 2X2 matrices 

Re(a^„-jS^„) 
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Mix 


Mil 
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A^31 


M32 


M33 





lm{amn+l^mn) 
- lm{(Xmn - ^mn ) Re(a«„ + ) 



(3) 



Here Re and Im denote the real and imaginary parts, respectively. The covariance 
matrix after a Bogoliubov transformation is given by = SgS^ . Let us suppose that we 
are only interested in the state of two modes k and k' after the transformation. A great 
advantage of the covariance matrix formalism is that the trace operation over a mode 
is implemented simply by deleting the row and column corresponding to that mode. 
Gonsider that the initial state of the quantum field is a general Gaussian state for 
modes k and /c', and that all other modes are in their vacuum state. The covariance 
matrix for modes k and k' is given by 

fkk' 

ij/k' and (^fcfc' are 2X2 matrices. Then the transformed covariance matrix is 



f^kk' = 



(4) 



where i/^/ 
given by 



where 



Ckk 
Ckk 



Ckk' 
Ckk' 



(5) 



Qj = Ml,iJy,Mkj+Mlf,,Mkj+Mlf,,^^^^ J] M^Mnr (6) 

Metrology techniques 

Having written the state of the field and its transformations in the covariance matrix 
formalism we are now ready to apply metrology techniques'^ that have been 
developed for continuous variable systems. In this section we briefly review some 
basic tools'^"^\ The aim in quantum metrology is to provide a strategy to determine 
the value of a parameter 0 with high precision when the parameter is not an 
observable of the system. Temperature, time, acceleration, and coupling strengths are 
good examples. A strategy corresponds to finding optimal initial states and mea- 
surements on the final states. In order to estimate the parameter with high precision it 
is necessary to distinguish two states p© and p@+d& that differ by an infinitesimal 
change d& of the parameter 0. The operational measure that quantifies the distin- 
guishability of these two states is the Fisher information'^ Let us suppose that an 
experimenter performs N independent measurements to obtain an unbiased 
estimator 0 for the parameter 0. The Fisher Information F(0) gives a lower 
bound to the mean-square error via the classical Gramer-Rao inequality^', i.e.. 



SCIENTIFIC REPORT: | 4 : 4996 | DOI: 1 0.1 038/srep04996 



2 



the likelihood function with respect to a chosen positive operator-valued measure- 
ment (POVM) {Ox} with 6a = • . Optimizing over all the possible quantum 
measurements provides an even stronger lower bound^^, i.e., 



iV( A© 



1 1 



(7) 



where H(0) is the quantum Fisher information (QFI). This quantity is obtained by 
determining the eigenstates of the symmetric logarithmic derivative Ap© defined by 



= ApQpQ + p^ApQ. Alternatively, the QFI can be related to the Uhlmann 



fidelity of the two states p© and +^^0 through 



H{0)- 



(8) 



where T{pi ,^2) = VfhPiVfh^ ■ The optimal POVMs for which the quantum 
Cramer-Rao bound (7) becomes asymptotically tight can be computed from A^©^^. 
Unfortunately, these optimal measurements are usually not easily implementable in 
the laboratory. Nevertheless, in typical problems involving optimal implementations 
one can devise suboptimal strategies involving feasible measurements such as 
homodyne or heterodyne detection, see, e.g.,^"^. Here we are interested in assessing 
metrology strategies based on the quantum Cramer-Rao bound and thus on the QFI, 
since our aim is to investigate how well one can in principle determine a parameter 
that appears in quantum field theory. As an example we will consider acceleration and 
we will show how our technique can be applied to develop a quantum accelerometer 
that exploits relativistic effects. By doing so, it is in principle possible to improve the 
state-of-the-art in accelerometer s. Before presenting our example we return to our 
general discussion. 

We consider a bosonic quantum field which undergoes a ©-dependent Bogoliubov 
transformation, where 0 is the parameter we want to estimate. For example, the 
transformation could be the expansion of the universe and the parameter in this case 
is the expansion rate. We assume that the initial state of the field is given by Eq. (4). To 
estimate 0 we must calculate the fidelity T{okh: [Q],Gkh: [0 + d&]), where the trans- 
formed state G]^k' [0] is given by Eq. (5). If g^}^ is a two-mode Gaussian state with zero 
initial first moments, the fidelity is given by^^ 



:F{dkkM,dkk\e+d^)-- 



n[0,0 + de] + ^n[0,0 + def - a[0,0 + ^^0] 



where n[0,0 + c/0] = v^r[0,0 + ci0] + v^A[0]A[0hM0], and 
r[0,0 + ^i0] = ^ det {Qokk' M^Okk [0 + ci0] - • ^ 



A[0] = -det(afcfc'[0] + /Q), 

A[0,0 + ci0] = det(^fcfc' [0] + ~Gk^ [0 + de]), 
io 



where • is the identity matrix. Note that we follow the conventions used in^3'24 ^j^^ 
normalization of the covariance matrix, which differ from other conventions^^. 
We now present an application of the techniques to measure accelerations. 

Results: Measuring Acceleration 

We consider a bosonic quantum field in flat spacetime confined in a 
cavity undergoing non-uniform motion. The motion-induced trans- 
formation of the quantum field depends on the magnitude of the 
acceleration. We will show that, using the techniques explained 
above, this fact can be used as the working principle of a highly 
accurate and sensitive accelerometer. 

Our aim is to determine the precision with which the acceleration 
can be estimated from measurements on the field modes. We present 
our analysis in (1 + 1) -dimensional spacetime with metric signature 
(— + ). Additional transverse dimensions can be included via their 
contribution to the effective field mass as discussed in^^'^^. We con- 
sider Dirichlet boundary conditions at the cavity walls. The details of 
the chosen boundary condition slightly modify the quantitative fea- 
tures of the modeP^, but not qualitatively. Physical implementations 
of this setup can be realised by optical cavities^^, superconducting 
circuits^^, or Bose-Einstein condensates^". 

The cavity is considered to be initially at rest and the state of the 
field given by Eq. (4). After a general trajectory the field modes 
undergo a Bogoliubov transformation. In order to treat the problem 



analytically we assume that the Bogoliubov coefficients which relate 
the initial and final state of the field have a series expansion in terms 
of a dimensionless parameter /z, such that 



where 



^ aL aLn^ 



(13) 
(14) 

(15) 



Here a is the proper acceleration of the cavity at its center, L is the 
length of the cavity in its instantaneous rest frame, is the propaga- 
tion speed of the excitations of the quantum field inside the cavity 
and we introduce n = cICs, where c is the speed of light in vacuum. In 
general the motion of the cavity can be an arbitrary combination of 
segments of uniform acceleration and inertial motion. For example, 
we can consider a finite period of uniform acceleration, a repetition 
of identical trajectory segments, or even sinusoidal oscillation with a 
fixed amplitude^^. 

The coefficients oc^^l = GmSmm where |G^| = 1, are the phases 
accumulated during both uniform acceleration and inertial seg- 
ments. Furthermore, the first order coefficients are zero on the diag- 
onal, i.e., a^^^ = P^^^ = 0. In this paper we suppose that the two modes 
k and k' have opposite parity, {k — k') is odd. This extra assumption 
causes the coefficients dy^^^' and [iy^h> to contain only odd powers of h in 
their series expansions. 

Let us consider a particular initial state of the modes k and k' . We 
assume the state to be two single-mode squeezed states in a product 

/ e^^'fc 0 

form. The 2X2 matrices in (4) are then i^^^ ( n ^-^r^ 



A9) ^,,= 



(10) 

(11) 

(12) 



0 



and (bj^y = 



0 e 

0 0^ 

0 e-^"^' J \0 0^ 

Currently, we are also investigating other initial states such as 
entangled Gaussian and Fock states. In this paper we restrict our 
analysis to the state mentioned above which already produces pos- 
itive results. The transformed covariance matrix Gkk'{h) is obtained 
using Eqs. (3), (5), (6) and (13), with our particular initial state. The 
QFI is given by 

8[1- 



^jT{~aaih),dkkih + dh))\ 



(16) 



where Tiokk' {h),dkk' {h + dh)) can be computed using Eq. (9). Using 
the Bloch-Messiah reduction [38, p. 9] it is possible to show that 
A = r + 0{h\dh^) eindA = 0{h\h + dhf). Therefore, the fidelity 

is givenhYJ^{dkk'{h),dkk'{h-\-dh)) = 1 — ^fti/z^, where f is the term 

proportional to dh^ in the Taylor expansion of F. Assuming for the 
sake of simplicity that both modes have the same squeezing para- 
meter Tk = Tk' = r, we obtain the QFI H{h) = H^^^ + H^^^ h^, where 



=5ft [4coshr(/^ +//) + 



4 cosh r a, 



- 4 cosher / 



4sinhV G-a>;! +Glp]^ +/; +/; _ Uu ; + 



2sinh2r| 24i>/;<i! +24Xi-cos(^*) I -/,'+/;-^ + 



) I -fi+fh 4^ + ^ ) +4sinhr(G-a:f + G-e-'',) 



[17) 



-4sinh2rcosh r 



-2sinh r a. 



sinh^2r 



r*2o/Vi 
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where 



0) 



and 



_^ ^ ^, a^^'* Pnj^*' The particular form of the Bogoliubov coefficients 

depends on the trajectory followed by the cavity. Arbitrary traject- 
ories composed of discrete intervals of accelerated and inertial 
motion were considered, for instance in^^, while continuous motion 
was addressed in^^. The latter case — continuous, sinusoidal motion 
with small amplitude — features two different kinds of resonances. If 
the frequency of the oscillation matches the sum of the frequencies of 

the two oddly separated modes the corresponding | | grows lin- 
early with the duration of the oscillation. This gives rise to a resonant 
particle creation phenomenon known as the dynamical Casimir 
effect^^. If, on the other hand, the oscillation frequency is equal to 
the difference between the frequencies of two oddly separated modes 

the resonance is associated to the coefficient |a[^)| (see ^^). In this 

paper, we take advantage of these resonances to increase the QFI in 
Eq. (17). Notice that our perturbative approach is restricted by the 
condition H^^^h^^l. 

Finally, using Eq. (7) the optimal bound to the error in the estima- 
tion of the parameter h after N measurements is obtained as 

(Ah) > ^ Assuming good control over the parameters L 

and c^, the final error in the estimation of the acceleration is just 
re-scaled by a factor of c] /L = c^ j (n^L) , that is: 

^2 



(Aa>> 



(18) 



Now we will consider a specific experimental implementation using a 
quasi one- dimensional BEC^^ on a flat spacetime metric^^'^^ with 
hard-wall boundary conditions^^'*^"^^. In the dilute regime, the EEC 
can be described by a mean field density plus phase fluctuations ^, 
which we expand in terms of the so-called Bogoliubov modes. The 
modes with frequencies well below the frequencies associated to the 
healing length of the condensate obey a massless Klein -Gordon 
equation in an effective curved spacetime metric. The effective metric 
depends on the real spacetime metric, the background pressure p, 
energy density p, number density n and flow velocity v^°'^\ In the 
absence of background flows v = 0 and for constant density (i.e. the 
homogenous case), the effective metric is also flat^^ 

Let us now describe the BEG after it undergoes acceleration, a 
situation that has been previously considered^^. In the comoving 
frame, the effective metric remains flat as long as the accelerations 
are small enough. In this case, we can ensure that the classical back- 
ground is not excited and that the condensate remains approximately 
homogenous. In^^ it has been shown that squeezing of the Bogoliubov 
modes occurs when a BEG with hard-wall boundary conditions 
undergoes small accelerations. The density of the BEG can in prin- 
ciple become inhomogeneous, however, these effects are negligible in 
the regimes considered in our discussion. Under the circumstances 
mentioned above, we obtain a Klein-Gordon equation for the 
Bogoliubov modes of the BEG with an effective flat metric and pro- 
pagation speed cl = c^dp / dp. The hard- wall boundary conditions 

fl c 

give rise to a mode spectrum given by C0n = 2nx , where L is 

the length of the condensate. Therefore, the techniques presented in 
the previous sections are directly applicable. A more detailed pre- 
sentation of our BEG setup can be found in^°. 

We have chosen to use a BEG to take advantage of amplification 
effects due to small propagation speeds among other convenient 
experimental features that we will discuss below. At this point, a 
comment on the relativistic nature of phonon production on a 
BEG is in order. In the Newtonian limit where v, Cg^c one can 
assume that Cs is independent of c as explained in more detail in 



reference^°. In this limit an accelerated observer and an inertial obser- 
ver are no longer related through Rindler transformations (i.e. 
Lorentz transformations in the uniformly accelerated case) but 
instead by Galilean transformations. In this case the time coordinate 
in the accelerated frame coincides with the time coordinate in the 
inertial frame since time dilation is negligible. As a consequence the 
cavity length remains constant in all frames and the vacua of the 
cavity at rest and the cavity undergoing uniform acceleration coin- 
cide. There is no particle creation, confirming that the effect we 
discuss is purely relativistic. In the case that v/c and cjc are not 
negligible, it is also possible to determine whether an effect is relat- 
ivistic or not. In this case an effect is considered to be relativistic if it 
disappears in the limit oo. By this criterion, photon production 
due to the motion of a boundary is a relativistic effect. In the BEG the 
excitations propagate much slower. However, the speed of propaga- 
tion is directly proportional to the speed of light through 
c] = dp / d p"^^'"^^ . Therefore, the effect also disappears in the limit 

^ 00 and we therefore consider it to be relativistic. 

We consider that the BEG undergoes sinusoidal acceleration given 
by /z(t) = h sin((X) t) and we use our scheme which exploits relativistic 
particle creation to estimate the amplitude h. For this continuous 
motion the Bogoliubov coefficients are given by^^. 



(t) = ie-' ^^^'^4'^) {(Dk' - (Ok) I dt sm{(D t)/^'^' -«^)^ 



4 



Jo 



(19) 



where aj.^) = 



3 and^[^) = - 



u^k'-kY ■ - + -tl^^ Bogoliubov 

coefficients that relate solutions to the Klein Gordon equation in the 
inertial and accelerated frames. It is easy to show that these coeffi- 
cients reduce to the identity when c^^co and in the Newtonian limit 
mentioned above. We assume that the frequency of the cavity oscil- 
lation is CO = C0k + (Dk' which generates a particle creation resonance. 
In Fig. (2) we plot the error in the estimation of the sinusoidal 
amplitude h for typical experimental parameters. The QFI, given 
by Eq. (17), for an initial squeezing of r = 10, mode frequencies cok 
= 2n • 500 Hz, w'k = 2wk, length L = 1 fim and effective velocity Cs 
= 10"^ m/s is approximately H(/z) ^H(a) ^ 10^^. Note that our 
perturbative analysis is restricted to a^lO~^m/s^. 

It is interesting to compare our methods with techniques that have 
been previously developed to estimate accelerations using BEGs 
within the framework of non-relativistic quantum mechanics. In 
particular we consider the QFI in state-of-the-art accelerometers 
designed for the detection of gravity anomalies^^'^^, which are based 
on Mach-Zehnder atomic interferometry. In these schemes the wave 
function of the BEG is split and recombined using laser pulses, giving 
rise to a phase shift ofcj) = k a T^, where k is the wave number of the 
atomic hyperfine transition, a is the average acceleration and T the 
interrogation time between pulses. The QFI in this case is given by^^ 
H = (da <t>y = (k V-y. Note that k is fixed by the atomic species 
employed {k = 1.6 • 10^ in Rb) and Tis limited by the dimen- 
sions of the experimental setup, being 1 s in the best scenarios. 
Therefore, H = 2.6 ' 10^^ m/s^. The readout of the phase is then 
performed by fluorescence imaging of the atoms. The optimal sens- 
itivity in the measurement of acceleration as provided by the QFI is 

given by^^'^^: ba=\ j (^VNkT^^ . After one measurement cycle, N is 

given by the number of detected atoms. After several cycles, the 
number of atoms is multiplied by the repetition rate and the integ- 
ration time. Gonsidering^^'^^ a number of atoms of 10^ — which 
already assumes a signal-to-noise ratio with respect to the total num- 
ber of atoms in the BEG — a repetition rate of 5 Hz and integration 
time of a few hours, a number N = 10^ Ms obtained. Therefore, those 
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Figure 2 | Optimal bound A a(m/s^) on the error in the estimation of the acceleration vs proper time of acceleration t{s) for squeezing parameter r = 8 
(blue, solid), 9 (red, dashed) and 10 (green, dotted). The frequencies of the modes are given by Wn = In ' 500 n Hz, as corresponds to L = 1 fim and 
Cs = 10"^ m/s. The motion is assumed to be a sinusoidal oscillation of maximum acceleration a and frequency co = cok cok' and the number of 
measurements is N = 10". The mode numbers are k = 1 and k' = 2. 



schemes require a stable large-atom BEC machine'*^ in order to 
achieve large integration times. The absolute sensitivity in^^'^^ is 
impressively small (Aa^5-10~^^m/s^). In comparable conditions, 
our scheme which has been designed within relativistic quantum 
field theory is capable of improving the QFI by two orders of mag- 
nitude. Assuming that it is possible to realize the same number of 
measurements as done in the non-relativistic schemes, the optimal 
bound to the absolute sensitivity in our setup is around 1*10"^^ m/s^, 
yielding a performance that goes several orders of magnitude beyond 
commercial accelerometers'*^ (Aa^ 10~^m/s^). 

The performance of our scheme can be further improved by 
employing entangled two-mode initial states or by increasing the 
number of measurements. As can be seen in Fig. (2), in our proposed 
implementation the error diminishes by increasing the degree of 
single-mode squeezing in the initial state. Notice that squeezing 
parameters of r = 10 appear to be within reach for phonons in 
time-dependent potential traps^°. Furthermore, a large number of 
measurements can be in principle achieved by using atomic quantum 
dots or optical lattices coupled to the condensate in order to probe 
the state of the phononic field^^. Assuming as in^^ that a few thou- 
sands of impurities can be coupled to the condensate and that each 
measurement can take a few milhseconds, the number of 10^^ can be 
achieved after several hours of integration -that is, in conditions 
comparable to^^'^^. We note here that even with a much more modest 
level of squeezing such as r = 2 and a number of measurements of 
10^, the predicted sensitivity of our device is 10~^, still two orders of 
magnitude beyond commercial accelerometers. Notice also that our 
setup is not restricted to a particular frequency of vibration, since 
there are resonances between any pair of oddly separated modes, 
both for particle creation and mode mixing. Therefore, a single 
BEC is sensitive to several frequencies. Moreover, each BEC can be 
tuned at will by changing the length of the trap or the speed of the 
propagation of the phonons, using standard experimental tech- 
niques^\ Taking all the above into account our setup can in principle 
exhibit a good broadband performance. 

Discussion 

The main aim of our research programme is the study of relativistic 
effects on quantum technologies. A comprehensive understanding of 
such phenomena will enable us not only to make the necessary cor- 
rections to technologies that are affected by them but also to use 
relativistic effects as resources. Indeed, we have shown that relat- 
ivistic effects can be exploited to improve quantum precision mea- 



surements. In particular, we showed how particle creation within a 
moving cavity, a quantum field theoretical effect known as the 
dynamical Casimir effect, can be used to determine accelerations 
with a precision that, in principle, can improve state-of-the-art in 
accelerometers. As a particular experimental implementation, we 
have discussed a BEC setup. We showed that the QFI is several orders 
of magnitude larger than its counterpart in non-relativistic schemes. 
Therefore, the ultimate bound to the sensitivity of the accelerometer 
is several orders of magnitude smaller. This means that, by employ- 
ing an optimal or close to optimal strategy of local phase estimation, 
our relativistic scheme can measure accelerations much smaller than 
the ones attained with optimal strategies in non-relativistic schemes. 
Moreover, we have presented a general framework that can be used to 
measure parameters that appear in quantum field theory such as 
gravitational field strengths, proper times and accelerations. 
Although Earth-based applications are also possible, the techniques 
are especially useful in space-based quantum technologies, where 
relativistic effects become relevant. 

This paper establishes relativistic effects as resources in quantum 
technologies. Our work opens an avenue for the development of a 
new generation of relativistic quantum technologies. 
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